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In the framework of the non-relativistic Brueckner-Bethe-Goldstone theory, we de- 
rive a microscopic equation of state for asymmetric and /3-stable matter containing 
E~ and A hyperons. We mainly study the effects of three-body forces (TBF's) 
among nucleons on the hyperon formation and the equation of state (EoS). We 
find that, when TBF's are included, the stellar core is almost equally populated by 
nucleons and hyperons. The resulting EoS, which turns out to be extremely soft, 
has been used in order to calculate the static structure of neutron stars. We obtain 
a value of the maximum mass of 1.26 solar masses (1 solar mass M ~ 1.99- 10 33 cj). 
Stellar rotations increase this value by about 12%. 

1 Neutron stars within the BBG approach 

The nuclear matter equation of state (EoS) is the fundamental input for build- 
ing models of neutron stars. These compact objects, among the densest in the 
universe, are indeed characterized by values of the density which span from 
the iron density at the surface up to eight-ten times normal nuclear matter 
density in the core. Therefore a detailecLknowledge of the equation of state 
over a wide range of densities is requiredu. This is a very hard task from the 
theoretical point of view. In fact, whereas at densities close to the saturation 
value the matter consists mainly of nucleons and leptons, at higher densities 
several species of particles may appear due to the fast rise of the nucleon chem- 
ical potentials. In our work we perform microscopic calculations of the nuclear 
matter EoS containing fractions of A andE - hyperons in the framework of the 
Brueckner-Hartree-Fock (BHF) scheme □. The BHF approximation, with the 
continuous choice for the single particle potential, reproduces closely the many- 
body calculations up to the three hole-line level. In this approach, the basic 
input is the two-body interaction. We chose the Paris and the Argonne Vis 
potential for the nucleon-nucleon (NN) part, whereas the Nijmegen soft-core 
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model has been adopted for the nucleon-hyperon (NY) potential. No hyperon- 
hyperon interaction is taken into account, since no robust experimental data 
are available yet. For more details, the reader is referred to ref. u and references 
therein. However, as commonly known, all many-body methods fail to repro- 
duce the empirical nuclear matter saturation point po = 0.17 fm~ 3 . This 
drawback is commonly corrected by introducing three-body forces (TBF's) 
among nucleons. In our approach we have included a contribution containing 
a long range_two-pion exchange attractive part and an intermediate range re- 
pulsive part □. This allows the correct reproduction of the saturation point. 
In figure 1 we show the chemical composition of /3-stable and asymmetric nu- 
clear matter containing hyperons (panel (a)) and the corresponding equation of 
state (panel (b)). The shown calculations have been performed using the Paris 
potential. We observe that hyperon formation starts at densities p ~ 2 — 3 
times normal nuclear matter density. The S _ baryon appears earlier than the 
A, in spite of its larger mass, because of the negative charge. The appearance 
of strange particles has two main consequences, i) an almost equal percent- 
age of nucleons and hyperons arc present in the stellar core at high densities 
and ii) a strong deleptonization of matter, since it is energetically convenient 
to maintain charge neutrality through hyperon formation than /3-decay. The 
equation of state is displayed in panel (b). The dotted line represents the case 
when only nucleons and leptons are present in stellar matter, whereas the solid 
line shows the case when hyperons are included as well. In the latter case the 
equation of state gets very soft, since the kinetic energy of the already present 
baryonic species is converted into masses of the new particles, thus lowering 
the total pressure. This fact has relevant consequences for the structure of the 
neutron stars. 



2 Equilibrium configurations of neutron stars 



We assume that a star is a spherically symmetric distribution of mass in hy- 
drostatic equilibrium. The equilibrium configurations are obtained by solving 
the Tolman-Oppenheimer-Volkoff (TOV) equations 111 for the pressure P and 
the enclosed mass m, 
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being G the gravitational constant (we assume c = 1). Starting with a central 
mass density p(r — 0) = p c , we integrate out until the pressure on the surface 
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Figure 1: In panel (a) we display the equilibrium composition of asymmetric and /3-stable 
nuclear matter containing £~ and A hyperons. In panel (b) the solid(dottcd) line represents 
the EoS obtained in the case when nucleons plus hyperons (only nucleons) are present. 
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Radius R (km) Central density n c (fm ) 

Figure 2: In panel (a) the mass-radius relation is shown in the case of beta-stable matter 
with hyperons (solid line) and without hyperons (dashed line). The thick line represents the 
measured value of the pulsar PSR1913+16 mass. In panel (b) the mass is displayed vs. the 
central density. The dotted line represents the equilibrium configurations of neutron stars 
containing nucleons plus hyperons and rotating at the Kepler frequency Q^- 
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equals the one corresponding to the density of iron. This gives the stellar 
radius R and the gravitational mass is then 

f R 

M G = m(i?)=47r / drr 2 p(r). (3) 
Jo 

For the outer part of the neutron star we have used the equations of state 
by Feynman-Metropolis- Teller □ and Baym-Pethick-Sutherlandtl, and ipr the 
medium-density regime we use the results of Negele and Vautherin 0. For 
density p > 0.08 fm -3 we use the microscopic equations of state obtained in 
the BHF approximation described above. For comparison, we also perform 
calculations of neutron star structure for the case of asymmetric and /3-stable 
nucleonic matter. The results are plotted in Fig. 2. We display the gravitational 
mass Mq (in units of the solar mass M ) as a function of the radius R (panel 
(a)) and central baryon density n c (panel (b)). We note that the inclusion of 
hyperons lowers the value of the maximum mass from about 2.1 M down to 
1.26 M a . This value lies below the value of the best observed pulsar mass, 
PSR1916+13, which amounts to 1.44 solar masses. However the observational 
data can be fitted if rotations are included, see dotted line in panel (b). In this 
case only equilibrium configurations rotating at the Kepler frequency Six are 
shown. 

In conclusion, the main finding of our work is the surprisingly low value of the 
maximum mass of a neutron star, which hardly comprises the observational 
data. This fact indicates how sensitive the properties of the neutron stars are 
to the details of the interaction. In particular our result calls for the need of 
including realistic hyperon-hyperon interactions. 
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